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^ ■ Abstract 
> 

, We study tlie metal-insulator transition in two sets of amorphous Sii_a;Ni^ 

I films. Tlie sets were prepared by different, electron-beam-evaporation-based 

technologies: evaporation of the alloy, and gradient deposition from separate 
Ni and Si crucibles. The characterization included electron and scanning 
tunneling microscopy, glow discharge optical emission spectroscopy, energy 
dispersive X-ray analysis, and Rutherford back scattering. Investigating the 
logarithmic temperature derivative of the conductivity, w = dlncr/dlnT, we 
I observe that, for insulating samples, w{T) shows a minimum, increasing at 

both low and high T. Both the minimum value of w and the corresponding 
O ■ temperature seem to tend to zero as the transition is approached. The analysis 

of this feature of w{T, x) leads to the conclusion that the transition in Sii_a;Nia; 
^ , is very likely discontinuous at zero temperature in agreement with Mott's 

?H ' original views. 
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I. INTRODUCTION 



Metal-insulator transitions (MIT) in disordered systems have attracted much interest, 
theoretical as well as experimental, for forty years. Milestones on this way were the concepts 
of Anderson localization of minimum metallic conductivity 0, the scaling theory of 
localization , and the renormalization group approach incorporating the electron-electron 
interaction into localization theory for surveys see Refs. 0^12 . 



Nevertheless, important questions have not been finally settled yet: Contrary to Mott's 
prediction of a finite minimum metallic conductivity, the scaling theory of localization pre- 
dicts continuity of the conductivity at the MIT for three-dimensional systems. That means, 
for the transition being approached from the metallic side at zero temperature, T = 0, the 
conductivity, a(T,x), is expected to vanish continuously according to 

cr(0, x) oc |a; — xd'^ , (1) 

where x stands for the control parameter (concentration / magnetic field / stress ...), and 
Xc for its value at the MIT. The value of the critical exponent u has been a subject of 
controversial debate since its first measurements, see e.g. Refs. |jl3|-|l8|. However, several 
studies have pointed to inconsistencies of the zero-temperature extrapolations inherent in 



the critical-exponent determination pJT3,hq-22|. These problems could, of course, arise from 



quantitative failure of the extrapolation formulae used or from imperfections of samples or a 
and T measurements, but they could also be caused by discontinuity of cr(0, x) at Xc [p|,p^, in 
which case the critical-exponent puzzle would simply be an artifact arising from unphysical 
fitting. The unclear situation described is the main motivation of the present study. 

The recent discovery of a sharp MIT in two-dimensional disordered systems [^3H 26|] has 
drawn interest to this subject, too. The scaling theory of localization of non-interacting elec- 
trons 0] states that two-dimensional systems exhibit principally only activated conduction, 
and denies the existence of an MIT in this situation. Hence it does not describe nature in 
this case. This failure is a strong additional motivation to re-consider the applicability of 
that theory to the MIT of three-dimensional systems. 

In greater detail, the critical step of most of the recent studies of the MIT in three- 
dimensional disordered systems is the zero-temperature extrapolation, that is the determi- 
nation of the limit a{T 0,x), for samples close to the MIT. This extrapolation can be 
based either on a microscopic theory, or on an empirically found relation supposed to be 
valid down to T = 0. 

The extrapolation is particularly difficult if a strong, but non-exponential T dependence 
of a is observed. In such a case, two interpretations are possible: 

(i) The sample could be metallic. This situation is mostly analyzed in terms of theories by 
Altshuler and Aronov P^, or Newson and Pepper [EH], which yield 



a(T, x) = a{x) + b{x) ■ (2) 

with p = 1/2, and 1/3, respectively. The former theory models the superposition of electron- 
electron interaction and disorder, but it is a perturbation theory so that its applicability very 
close to the transition is at least questionable. The latter theory considers the T-dependent 
drop of the diffusion constant as the decisive variation, and yields a power law with exponent 
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1/3 for X = Xc- 

(ii) The sample could exhibit activated conduction with the characteristic energy being 
smaller than the lowest experimentally accessible T. Thus it would have to be classified as 
insulating (at T = 0). If, as generally taken for granted, the characteristic energy vanishes 
continuously as a; — > Xc, such an x region exists always, see Fig. 1. But, to the best of our 
knowledge, there is no appropriate microscopic theory for a related quantitative analysis 
available. 

Unfortunately, when considering non-exponential cr{T, x = const.) dependences, many of 
the recent experimental studies analyze the data only in terms of (i). Moreover, to agree with 
Eq. (1), they presume that, as x — > Xc, the parameter a tends to 0. So it is not surprising 
that, when checking the applicability of Eq. (2) by considering the logarithmic derivative 
dlncr/dlnT, doubts arose in several cases |T9|-p^]: this quantity must vanish as T — >• 
for metallic samples, but, instead, often increases or stays approximately constant for the 
sample(s) assumed to be particularly close to the MIT. Therefore one should determine the 
MIT point by simultaneously considering both sides of the transition. Moreover, adjusting 
free parameters should be avoided to the largest possible extent. Such an approach was 
successful for amorphous Sii-a^Cra,, see Refs. p9|-|33|. A simple phenomenological model 
of o"(T, x) describing both sides of the MIT was constructed for that system starting from 
the detection of universal features in a{T,x), in particular from a scaling law for the T 
dependence in the activated region, a(T,x) = a{T /Tq{x)). However, annealing breaks this 
universality. 

Due to the controversial situation described above, the aim of the present work is to 
analyze the (t(T, x = const.) relations close to the transition in as unbiased a way as possible, 
without fitting. For that we study amorphous silicon nickel films, a-Sii_2;Nij., prepared by 
electron beam evaporation. Thus, though the above discussion concerns a general problem 
occuring for an arbitrary type of the control parameter x, we focus now on a special case, 
and, in the following, x stands only for metal concentration in an amorphous alloy. In our 
study, we put emphasis on two methodical points: Firstly, we attach particular importance 
to the accuracy of the dlncr/dlnT values. Secondly, in order to avoid being mislead by 
having used special preparation conditions, we investigate two series of samples, prepared 
by groups in Dresden and Tel Aviv. 

There were several reasons for us to select a-Sii-^^Ni^^ for a detailed investigation: 
- We were looking for an alloy which is sufficiently stable but differs substantially from 
a-Sii_a;Cra, which we had studied in detail previously [P^-|55[]. According to the melting 



temperatures ||3^, a-Sii-^^Nia; can be expected to be rather stable, at least far more than 



a-Sii_a;Au2; [^, or a-Gei^x-^Ux [0. On the other hand, there is a deep eutectic point 
close to the MIT for the system Si-Cr, but not for Si-Ni [Q. Therefore, in a quenching 
experiment, it should be far more difficult to realize an amorphous structure close to the 
MIT in Sii-i^Ni^: than in Sii_i.Cr^. |^. However, it is an open question, to what extent such 
considerations can be applied also to films produced by evaporation, or by sputtering. 
- A second important difference between a-Sii_i.Ni2. and a-Sii_a;Cr^. concerns the character 
of the silicide which should be formed first in annealing the samples, cf. Refs. p8|-^: Si2Ni 
is a metal, whereas Si2Cr is a semiconductor with an indirect gap of roughly 0.3 eV [ET . 



- Conductivity studies (partly taking magnetoresistance into account) on hydrogenated and 
unhydrogenated sputtered a-Sii_a;Nia, films as well as on polycrystalline films P^-HS 
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and on the influence of annealing [|^,^ were available for comparison in the literature. 
Moreover, in the vicinity of the MIT, the electronic structure of a-Sii-a^Nij. was recently 
studied by XPS and UPS 0]. For a-Sii_^Nia.:H, a structure study |511], optical conductivity 



measurements also under pressure |5^, and electronic structure investigations by XPS 



and XES were reported in the literature. Last but not least, some of the authors have 
gained detailed experience in trying fits of cr(T) by conventional theories for a-Sii-ajNi^; 
samples prepared by electron-beam evaporation in Tel Aviv | pB| , |5B |. 

The present paper is organized as follows: Section II describes the sample preparation. 
Section III is devoted to the structural and chemical characterization, including an STM 
investigation. Section IV gives some experimental details of the conductivity measurements, 
whereas Section V analyses them in detail, and lists common and differing features of the 
two conductivity data sets. In Section VI, we discuss the qualitative character of the MIT. 
Finally, in Section VII, the conclusions from our findings are summarized. 



II. FILM PREPARATION 

We compare a-Sii-^Ni^. films, prepared by two different technologies. Both technologies 
are based on electron-beam evaporation, but they have specific advantages and disadvan- 
tages: 

(A) Direct evaporation of SiNi alloys of various compositions from one crucible was used 
in Dresden, cf. Ref. [^]. In this case, changes of power cause changes of the composition 
of the films by modifying the vapor pressure ratio of the alloy components via the changed 
temperature of the ingot. Hence, the deposition rate has carefully to be held constant. 
Moreover, it is essential that the whole ingot is liquid; a check for contamination from the 
copper crucible showed it to be negligible. The x drift over the film thickness arising from 
Ni enrichment in the ingot material during evaporation is very small, see section III, since 
only a small portion of the ingot was evaporated in each cycle. Moreover, this drift is far 
smaller than the x difference between samples prepared in successive cycles using the same 
ingot, because the time which we waited until rate and pressure were stable is considerably 
larger than the deposition duration. Typical conditions were as follows: residual gas pres- 
sure during deposition ^ 4 ■ 10~^ mbar, and rate ~ 2 nm/s. 

(B) Co-evaporation of Ni and Si from separate crucibles was used in Tel Aviv to prepare in 
one deposition a whole set of samples. This method uses the varying incidence angles and 
distances to the sources to produce a composition gradient along a substrate. The film is 
then cut into narrow samples with their axes (the direction of current flow) perpendicular 
to the composition gradient This technique has two advantages: Composition deter- 



mination is related to a geometrical measurement, and neighboring samples are prepared 
under almost identical conditions. Thus the random errors in studying x dependences can 
considerably be reduced. However, the crucial point is that one has to take great care to 
hold the ratio of the evaporation rates of both crucibles constant. The conditions were as 
follows: residual gas pressure during deposition ^ 1 ■ 10~^ mbar, and, for the samples close 
to the MIT, rate ~ 0.8 nm/s; for more details see Ref. . 



In the following, we refer to the samples prepared by technology A by the numbers 1 - 
6, and to the samples produced by method B by the letters a - j. These two sample sets are 
labeled as A and B according to the preparation. Tables I and II give x, film thickness, and 
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conductivity at 300 and 4.2 K for the sets A and B, respectively. Samples 6 and i, both deep 
in the metallic region, and sample j, close to the MIT, are only included for comparison of 
X scales. 



III. CHEMICAL AND STRUCTURAL CHARACTERIZATION 

A. Composition 

The chemical composition of the samples of set A was determined by Rutherford 
backscattering (RBS) at the Research Center Rossendorf within this work. For that we 
used Sigradur substrates (glassy carbon). We confirmed the consistency of these data by 
glow discharge optical emission spectroscopy (GDOES). Energy dispersive X-ray analysis 
(EDX) by means of a LINK ANIOOOO EDS system attached to a JSM-840 SEM was used in 
determining the composition for sample set B in Witwatersrand, South Africa, within Refs. 
| 55| , |56[] . The results of these analyses are given in Tables I and II. 



In order to compare both x scales, we performed additional EDX analyses for samples 
3, 6, f, h, and i by means of a Tracor/Noran spectrometer Voyager Ila, attached to a 
transmission electron microscope (TEM) Philips CM20FEG, in Dresden. For that, pieces 
of the film were scraped from the substrate with a steel knife, and placed on a carbon-film- 
coated copper microscope grid. Moreover, for sample j, the Witwatersrand EDX result was 
checked by an RBS analysis performed in Faure, South Africa, within Ref. [ p^ . 

The results of these analyses are presented in Table III. It illustrates the inconsistencies 
in the x scales which were already pointed out in Ref. ||5^. Note that the Dresden EDX 



data underestimate our Rossendorf RBS values by roughly 5 at% Ni for sample 3 as well as 
for sample 6. Moreover, comparing the Dresden EDX data with the Witwatersrand EDX 
values we conlcude that the latter overestimate the Ni content, according to our RBS, 
by roughly 6 at% close to the MIT. This is confirmed by the comparison of South African 
EDX and RBS data for sample j, yielding a difference of 5.3 at%. For sample i, the EDX 



value of Ref. ||55| might be very close to the Ni content according to our RBS scale - Table 
III suggests that the systematic deviation of the Dresden EDX values from our RBS data 
could be roughly constant. Finally, an additional check with an SiNi2 bulk standard showed 
the Witwatersrand EDX to overestimate the true value by roughly 3 at% there. 

We consider the RBS as likely to give the most accurate x values since this method 



is quantitative from first principles and does not require elemental standards This is 
confirmed by the good agreement of both RBS scales in Table III. In our case, the relative 
uncertainty of the corresponding Ni:Si ratios was estimated to be better than 2 %, so that 
the related absolute error of the Ni content should not exceed 0.3 at% in the MIT region. 
The quantitative result of the EDX analysis however is influenced by a series of effects, in 
particular this method needs standards and various corrections . Therefore, the system- 
atic errors can amount to several at% so that the differences between the columns of Table 
III are not surprising. 

The X data, given in Tables I to III, and discussed in the previous paragraphs, are 
obtained from the Ni:Si ratio ignoring all kinds of contamination. Different kinds of analysis 
showed oxygen contamination to be the most relevant one. By means of EDX, we obtained 
4.4 and 11 at% as oxygen concentration in the samples 3 and h, respectively. If the oxygen 
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were distributed homogeneously in the films, a Ni concentration of 20 at% determined from 
the Ni:Si ratio would correspond to true 19.2 and 18 at% Ni for series A and B, respectively. 
However, some of the oxygen is expected to be concentrated close to the film surface and 
the film-substrate interface, so that this quantitative information should not be overvalued. 

Both sets of films were found to be well homogeneous. Laterally we checked this by EDX 
and electrical measurements for both sets. Changes in vertical direction might be critical 
in both technologies. We used GDOES to clarify this point. Because of the non-conducting 
substrate, we applied our self-developed radiofrequency equipment, a new and very suitable 
device for studying vertical concentration variations [^. For series B, the vertical homo- 
geneity was also checked by monitoring the quotient of the Ni and Si evaporation rates 
during film deposition [l55| . 

By means of GDOES, the Ni:Si ratio was found to deviate from the respective mean 
values through the thickness of the samples by not more than 5 % for both sample sets, 
corresponding to ±0.7 at% Ni close to the MIT. Fig. 2 shows depth profiles for series A and 
B. Note the constancy of the Ni signal. The variation of the Si signal at the beginning is 
an artifact. It arises from the strong dependence of the intensity of the Si line used on the 
changing plasma parameters at ignition. This variation is caused by surface contamination 
indicated by the C (carbon) signal; this was checked by additional experiments. Therefore, 
we estimated the relative variation of the Ni:Si ratio from the Ni signal. 



B. Structure 

The microstructure of the samples of both sets was investigated by TEM. In all cases, 
the samples were found to be well amorphous. However, the existence of crystalline regions 
with diameters below 2 nm cannot be excluded. Such a clustering is suggested by cf. 



The topography of the a-Sii-^Ni^; films was studied ex situ by means of a scanning tun- 
neling microscope (STM) under UHV conditions {p < 1 ■ 10~^mbar). For that purpose, we 
used a modified Omicron-UHV-STM with etched tungsten tip and partly self-made electron- 
ics. A typical topography of sample 3 is given in Fig. 3. Three corresponding cross-sections 
are presented in Fig. 4. For sample f, we obtained very similar results. 

Figs. 3 and 4 show that there are fluctuations in the film thickness of typically ±2 nm. 
In rare cases, the fluctuations reach an altitude of ±7 nm. The horizontal diameter of these 
hills amounts to the order of 30 nm. However, there is also a fine structure on a smaller 
scale, which fact might be related to some fractal structuring Unfortunately, we have 
no information to what extent these structures can be considered as representing a surface 



roughness only, or whether they originate from a columnar film growth, see e.g. Ref. ||6^ . 

The surface roughness could cause differences between technologies A and B concerning 
the chemical homogeneity of the films when considering a mesoscopic length scale: The local 
chemical composition is independent of the surface slope in technology A. However, during 
gradient film deposition by means of technology B, the variation of the surface slope causes 
a fluctuation of the Si and Ni incidence angles, where an increase of the former is correlated 
with a decrease of the latter. Thus the chemical composition of these films ought to fluctuate 
on a mesoscopic length scale |^ . 
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To estimate qualitatively this effect is difficult for two reasons: (i) The ffuctuation 
strength depends on the size of the surface area taken into account in averaging, (ii) We 
have no information on surface and volume diffusion lengths limiting composition ffuctua- 
tions. As an example, we consider a cube with edge length 2.4 nm, and assume the diffusion 
lengths not to exceed this value. The typical surface slope, estimated as root mean square 
of the difference quotient Az/Ax (with Ax = 2.4 nm) for several cross sections, corresponds 
to an angle of 14 degrees. For sample f, the incidence angles of the Si and Ni atoms, arriving 
from opposite sides, are roughly 10 and 45 degrees. This sample has a mean Ni content 
of 24.8 at% according to the EDX scale of Ref. Thus the ffuctuation of the incidence 
angle due to the surface roughness causes composition ffuctuations with a mean deviation 
of ±6 at% Ni. This value is by a factor of 3 larger than the finite size induced random dis- 
persion of the Ni content. However, since the regions of maximum and minimum Ni content 
do not form a percolative network, they have to be considered as local inclusions. 



IV. CONDUCTIVITY MEASUREMENTS 

Three types of measurements have been performed to study the T dependence of a: 

(a) 2 - 300 K using a self-built cryostat insert based on a bell technique, which permits very 
accurate stabilization of T over the whole interval, so that a high accuracy of the values of 
the logarithmic derivative is guaranteed, 

(b) 0.45 - 4.2 K using an RMC ^He cryostat, and 

(c) 35 - 900 mK using a dilution refrigerator Kelvinox 300 from Oxford Instruments. 
Measurements (a) of the B samples are repetitions with increased accuracy of the inves- 



tigations in Refs. |]55| , ^ . Some of the measurements (b) were already published in Refs. 
55| , |56[] . However, within this study they were re-analyzed, and the values of the logarithmic 
derivative were re-calculated by means of Eqs. (8,9) below. 

Two methods of making contacts to the samples were used: The leads were attached 
either by pressing indium tabs onto the sample, or by silver paste. Both these technologies 
avoid post-preparation heating of the films, but they need care to avoid contact problems 
in thermal cycling. 

A specific problem consists in the aging of the samples which can cause a resistivity 
change up to the order of 10 %. Thus we performed annealing experiments: Metallic films 
are infiuenced only slightly by annealing below roughly 150°C, but annealing up to about 



300°C leads to great changes, cf. Refs. P^JI5[ ]. Repeated thermal cycling between cryogenic 
and room temperatures can cause resistivity changes by a few % as well, but no drift is 
detectable in thermal cycling between 2 and 4.2 K. The usual aging concerns only the 
absolute value of the resistivity; the resistivity ratio R(T) / R{300 K) was observed to remain 
almost unchanged. Therefore we do not consider the sample after aging / thermal cycling 
to be a different sample, as if prepared under different conditions, and we simply scale the 
resistivity when necessary to retain the original value at an appropriate T. However, this 
scaling is irrelevant for the consideration of the logarithmic derivative, on which the main 
conclusion of this paper is based. 

Overviews for both the data sets A and B are shown in Fig. 5. Both data sets exhibit 
consistent behavior, with some systematic differences between them. For comparison. Fig. 
5a includes data of two sputtered a-Sii_^Ni^. films from Ref. |HD[, resembling the curves 
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of our evaporated samples very much. However, it is almost impossible to obtain reliable 
information on physical mechanisms directly from these only weakly structured curves. 



V. CONDUCTIVITY DATA ANALYSIS 

A. Standard extrapolation 

The standard analysis for determining the critical concentration of the MIT, based on 
extrapolation according to Eq. (2) with p = 1/2, is presented in Fig. 6. Indeed, if only data 
points between roughly 2 and 30 K are taken into consideration, linear relations correspond- 
ing to Eq. (2) are found in the a versus T^/^ representations. According to this analysis 
the MIT should be located just below the x of sample 1 in set A, and between samples c 
and d in set B. However, inclusion of the data around and below about 1 K shows that this 
classification is questionable. 



B. The logarithmic derivative 

More insight can be gained by considering the logarithmic derivative, 

(3, 

Its limiting behavior for vanishing T allows one to detect the transition point ||19 |. 



As a first approach to an understanding of how w{T, x) reflects the character of the MIT, 
let us assume that, for a given sample, the logarithmic derivative of the conductivity can be 
approximated by 

w{T)=c + d-{T/T,,,y , (4) 

where Tm is an experimentally accessible temperature, at which a has been measured. 
Integration of Eq. (4) yields for q ^ 0: 

In a(r^) - In a(T) = c ■ (In - InT) + d ■ [1 - {T/T^y]/q , (5) 

so that 

a{T) = a{T^) ■ {T/T^Y ■ exp{-rf[l - {T/T^Y]/q} . (6) 

For g = 0, we can assume d = 0, and get from Eq. (4): 

a(T) = a{T^) ■ {T/T^f . (7) 

Interpreting Eqs. (6,7), we have to discriminate between the eight situations listed in Ta- 
ble IV. According to this table, w(T — > 0) = indicates metallic character of the conduction, 
whereas both diverging and finite positive w(T ^ 0) point out that the sample is an insula- 



tor, cf. Refs. p!9|- p^ , |6^ . In more detail, diverging and finite positive w{T — > 0) are correlated 



with vanishing of a as T ^ according to exponential and power laws, respectively, cf. Eqs. 
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(6,7). Studying an MIT, we can exclude the case that the hmit w(T — > 0) < 0, correspond- 
ing to a(T — i> 0) = oo, as unphysical. However, this does of course not exclude negative w at 
finite T, which can arise from weak antilocalization, or e.g. from electron-phonon scattering 
in good metals. 

In this context, we may ask what the condition is for observing exponential behavior, 
(t(T) oc exp{ — (To/T)"^}. The condition T < Tq implies w > q. Thus there exists a low-w 
range where, though increase of w with decreasing T would rule out metallic conduction, it 
would still not be possible to identify activated conduction by detecting exponential criT). 
However, assuming only monotonicity of if as T ^ 0, one reaches an unambiguous decision. 

The logarithmic derivative is also helpful in identifying the physical mechanisms involved 
p9| , |30|j65| since it makes characteristic features of ct{T) visible, which are not obvious from 
(t(T) itself. An additional advantage of this quantity is that it is not infiuenced by geomet- 
rical errors, nor by the scaling of the resistivity mentioned above. 

Our method of data evaluation for analyses in terms of w is based on that used in Ref. 
|66| : A linear fit of In a versus In T for k neighboring points is performed. 



_ fc 'Ei 111 (^i 111 Tj - Y.i In (Jj E j 111 Tj 
fcE.(lnT.)2-(E,lnT,)2 " 

In order to minimize the numerical error arising from the non-linearity of the dependence 
of In cT from InT, we relate this w value to a mean temperature, Tif, defined by 



2(A:E.(lnT.)2-(E,lnT,: 



Thus w\{ would exactly equal w{Tif) if In a were a second order polynomial of In T. Utilizing 
Eqs. (8,9) allows one to take into consideration rather broad T intervals, so that the total 
error (numerical plus random) can be kept small. 

The optimum value of k depends on the accuracy and density of the data points, as 
well as on the structure of the measured curve. We considered the following groups of 
neighboring data points: The measurements in the region 2 - 300 K yielded data points 
with comparably high accuracy, but low density. From them, the logarithmic derivative 
was obtained using k = 3 below 7 K, and k = 2 above. In the ^He region, the density 
of data points is considerably higher, thus 8 neighboring data points were considered in 
differentiating for all samples except for h, for which we used 4 points. The data obtained 
in the dilution refrigerator were evaluated using k = 5. 

Fig. 7 gives a survey of the T dependences of the logarithmic derivative. Low-T details 
become visible in Fig. 8. These figures show that both the sample sets A and B are internally 
consistent. However, there are also significant differences between these sets. 



C. Common features of the two series 

Consider first the common features: For both data sets, at sufficiently high T, all samples 
studied behave very similarly: w{T) increases with increasing T. On the other hand, for the 
samples with the smallest Ni content, at low T, there is a pronounced increase of if (T) with 
decreasing T, indicating activated conduction. The strength of this contribution decreases 
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with increasing x. One common feature of all samples (with w > 0), which exhibit such 
an increase of w with decreasing T, has to be stressed: this increase continues down to 
the lowest accessible T. The extrapolation / generalization of this finding will be a basic 
assumption of our data analysis. 

The differing behavior at low and high T must be caused by two different mechanisms 
being dominant at low and high T, respectively. The low-T mechanism related to increasing 
w{T) as T decreases is very likely a kind of hopping conduction, cf. Refs. P^,pD|,|SB|,pB 



Concerning the origin of the high-T contribution we can only speculate. Such a high- 
T contribution is common to many amorphous transition- metal semiconductor alloys [Q. 
Comparing energy scales, it was speculated that in a-Sii-^^Cr^^ this mechanism might be 
related to electron-phonon interaction |^ . 

The feature, which we consider to be particularly interesting in Figs. 7 and 8, is the 
minimum of w{T) related to the crossover between these two mechanisms. This minimum 
is located at T^in = 150 K for sample a, which is the most insulating sample studied. The 
related minimum value of w is Wmin = 0.68. The minimum shifts to lower T with increasing 
Ni content. This behavior has to be expected since the characteristic hopping energy should 
tend to as the transition is approached, so that hopping only becomes dominant at lower 
and lower T. 

Simultaneously with T^in, "U^min decreases: For the A samples 1 and 2, T^i^ = 6 K and 
4 K, respectively, where Wmin = 0.42 and 0.32. For the B samples d and e, T^in = 0.8 K and 
0.2 K, respectively, where Wmin = 0.15 and 0.06. 

The four samples 1, 2, d, and e must be insulating because, for them, w(T) increases with 
decreasing T at the low-T end of the T range considered. This increase cannot be understood 
in terms of metallic behavior according to Eq. (2). Moreover, because of the overall behavior 
of w(T, x) (see first paragraph of this subsection), it is very likely that this trend continues as 
T —>■ what indicates insulating character according to Tab. IV. Furthermore, even if w{T) 
does not increase as T — > but only remains constant, i.e., w{T < Tmin) = Wmin, the samples 
still have to be classified as insulating, because a{T < T^in) oc T'"""", cf. previous subsection. 
This finding rules out the classification of these four samples as metallic according to the 
standard procedure used in Fig. 6. 

However, for the samples 3, 4, and 5, w{T) seems to tend to as T ^ 0, as well as for 
the samples g and h. (In this judgement, we presume that, since we study an MIT, the limit 
w(T ^ 0) > 0, see previous subsection.) For sample f, we cannot decide whether or not 
w{T ^ 0) = although we investigated it down to 35 mK. Thus we conclude that the MIT 
should take place between samples 2 and 3 for set A, and between samples e and g for set 
B. 

In the calculation of the critical Ni contents, and of their possible errors, we take into 
account the Ni concentrations of the samples bracketing the transition, the uncertainty of 
the X values, and the influence of preparation conditions differing slightly from sample to 
sample in technology A. For sample sets A and B, we obtain 16.6 ± 1 at% according to our 
RBS scale, and 25 ± 2 at% according to the EDX scale of Ref. respectively. In order 
to compare both values, we scale the EDX data by means of Table III: Presuming a linear 
relation between both EDX scales close to Xc, and the RBS value to exceed the EDX value 
of this study by 5 at%, one obtains 19 ± 3 at% as the critical Ni content of sample set B 
according to our RBS scale. Additionally performing the correction concerning the oxygen 
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concentration, see Section III, we obtain similar values for the critical concentration: 15.9±1 
and 17 ± 3 at% Ni for series A and B, respectively. Note however, that this agreement is 
not a necessity. On the contrary, since O is very reactive, it would not be surprising if the 
O contamination significantly influenced the gap states and shifted x^, cf. the role of H in 
a-Si based alloys p7| . 

The consideration of the samples bracketing the transition in both sets shows that, at 
high T, samples with x just below and above Xc exhibit very similar T dependences of w - 
there is no experimental indication of a discontinuity in the x dependence at finite T. Hence 
the respective substantial variations of a with T on both sides of the MIT must originate 
from the same physical mechanism, which, crudely speaking, "survives" the transition. This 
high-T mechanism yields a substantial contribution to w down to the lowest T studied, 
compare the samples closest to the MIT. Thus there is a superposition of the low- and high- 
T contributions, whatever mechanisms they represent, in the whole T interval investigated. 
This has an important consequence: In fitting simply to a pure hopping relation close to the 
MIT, or to Eq. (2), one might obtain a precise mathematical description, but the derived 
parameters would not have physical relevance. 

Concerning the high-T mechanism, the comparison with a-Sii-xCr-c given in Fig. 9 is 
interesting. For both the a-Sii_a,Cr^ sample 6 and the a-Sii_a;Nij,. sample 4, w(300 K) ^ 0.5. 
However, with decreasing T, the influence of the high-T mechanism decreases considerably 
faster in sample 6 than in sample 4: w{30 K) = 0.04 and 0.16 for samples 6 and 4, respec- 
tively. This faster decrease of the influence of the high-T mechanism is a general feature of 
a-Sii-xCri:. It has two consequences: (i) Comparing insulating samples of both substances 
with the same Wmin, the minimum is reached in aSii-xCvrc at clearly higher T than in a- 
Sii_j;Ni;j., compare sample 7 with sample 2 in Figs. 9 and 7a, respectively, (ii) In a-Sii-^jCrj;, 
there are minima in ct(T), i.e., zeros in w{T), at significantly higher T than in a-Sii-ajNi^;, 
consider samples 6 and e, and see next subsection. Thus, in contrast to a-Sii-^^Nia;, there 
is a broad, easily accessible low-T region in a-Sii-xCr^;, where one can ignore this high-T 
contribution, and study the pure low-T mechanisms. 



D. Differences between the two series 

Besides the common features described above, there are significant differences between 
the data sets A and B. These differences concern both the high-T and the low-T regions. 

Comparing A and B samples with similar low-temperature w{T) values, the A samples are 
found to exhibit larger high-temperature w values, i.e., they are more strongly T dependent. 
In detail, w(300 K,Xc) = 0.67 ± 0.06 and 0.48 ± 0.06 for sets A and B, respectively. This 
might indicate a similarity to a-Sii-xCr^;, where the high-T mechanism seems to lose its 



influence with increasing oxygen content of the films, compare Fig. 2 of Ref. |]68[ with Fig. 1 
of Ref. The same seems to be valid for a-Sii_a;Nia.: The high-T contribution has smaller 
influence for series B, for which the oxygen contamination is larger, as observed by EDX. 

For samples exhibiting a decrease of if (T) as T — > 0, this decrease is steeper for the A 
than for the B samples. Sample 3 even exhibits a zero of w{T) at about 0.9 K, and for 
sample 4 a zero just below the lowest measuring T is very likely. But such a zero was not 
found for the B samples within the T region studied. 
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Thus, in comparing the metalhc samples of both sets, i.e., 3, 4, 5 and g, h, it is puzzhng 
that for the A but not for the B samples there is visible a contribution to o"(T) with negative T 
derivative. Under what conditions such a contribution exists is a question of current interest 
|69| . E.g. such contributions were observed in a-Sii_a;Cr;j. [Q, as well as in heavily doped 



crystalline Si:P and Ge:Sb [^. It remains an open question whether the missing of 
this contribution in the B samples is real, or whether this contribution is obscured in set 
B by the low-T part of the high-T contribution. The latter hypothesis is suggested by the 
comparison in Fig. 7a which demonstrates that this low-T part is considerably larger in the 
B samples than in the A samples. 

Because our experimental data on the contribution with negative da/dT are very sparse, 
we can only speculate on the physical mechanism behind. It is likely that its origin is 
the same as in a-Sii_a;Cra: p0|,|32|, see also Fig. 9, heavily doped crystalline Si:P |70], and 



Ge:Sb [^, that means probably some interplay of electron-electron interaction and disorder. 



However, the interpretation in terms of precursors of a metal-superconductor transition as in 



a-Sii_a;Au2;, where both the pure constituents do not become superconducting ||3^, cannot 
be ruled out yet. Further investigations at lower T are desirable. 

Since great care was taken for both sample sets concerning homogeneity, we ascribe the 
differences stressed to the influence of differing preparation conditions. Significant contam- 
ination by crucible material could be excluded for both sample sets: In a related GDOES 
study, the contamination of an A sample by copper from the crucible was found to be below 
500 ppm. Concerning a crucible-material-contamination check for the B samples see Ref. 



| 55| . Thus only three possibilities seem to remain: 

- Contamination by any constituents from the residual gas, 

- differing surface diffusion during deposition arising from the differences in rate and / or 
temperature, and 

- concentration fluctuations on a mesoscopic length scale arising in the technology B, cf. 
Section III. 

Further investigations are needed to clarify this point. 



VI. DISCUSSION: THE QUALITATIVE CHARACTER OF THE MIT 

A. Missing plateau in w(T) 

The common features of both sample sets suggest an important qualitative conclusion 
concerning the character of the MIT, as will be explained in this section. Figs. 7 and 8, 
together with Tables I and II, show that w{T = const., x) decreases monotonically with 
increasing x, at least as long as w > 0. Thus there must be a separatrix (separating line), 
above which insulating and below which metallic samples are located in these graphs. More 
accurately, instead of a separatrix, there could also exist a separating strip, but this would 
imply a discontinuity in a{T = const., x) for which there is no experimental evidence in this 
experiment, nor in the literature. 

Now we compare the consequences of the common assumptions with experiment. For that 
we assume that the transition is continuous, as described by Eq. (1), and that, sufficiently 
close to the transition, below some temperature T*, Eq. (2) holds. These equations imply 
three conclusions on w{T,x): 
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(i) w{T,Xc) — p since a{xc) — 0. Thus the separatrix should be parallel to the T axis 
below T*. 

(ii) For metallic samples close to the MIT, w ^ p as long as a{x) <C b{x) ■ T^. In other 
words, since a{x) vanishes the logarithmic derivative w{T,x — const.) should 
exhibit a plateau below T*, which extends to lower T the closer x to Xc- 

(iii) A plateau should also be present on the insulating side close to Xc below T* for the 
following reasons. Since the charatcristic hopping energy very likely tends continuously to 
as X ^ Xc, the crossover temperature, below which the influence of hopping dominates the 
T dependence of w, vanishes too. Continuity of a{T — const., a;) leads to the expectation 
that w ~ const, above that crossover temperature. 

However, according to our experimental data, such a plateau is missing. Instead, for 
arbitrary fixed T, w{T,Xc) must be smaller than w{T) of samples 2 and c, but exceed w(T) 
of samples 3 and g for sample sets A and B, respectively. Thus, w{T,Xc) decreases with 
decreasing T, at least down to 0.2 K, where w{T,Xc) < 0.06. There is no experimental hint 
for w{T, Xc) at even lower T to rise again up to one of the theory values 1/2 or 1/3, and to 
saturate there; for x close to Xc, such a behavior would imply the existence of metallic tv{T) 
curves with a low-T maximum, so that it would contradict the assumption on monotonicity 
of w{T) as T vanishes for samples with w > and dw/dT < 0, concluded in Section V.C 
from our experimental data. Hence, at least one of the standard descriptions Eqs. (1,2) is 
very likely not valid: Either the relevant mechanism at the transition is related to a value 
of p which is far smaller than expected, or the transition is not continuous at T = 0, since 
w{T 0, Xc) — indicates finite a{T 0,Xc + 0) according to Table IV. 



B. Two simple models 

To illustrate the problem from a different perspective, we study w{T,x) for two qualita- 
tive empirical models. They are constructed so that continuity at Xc for arbitrary finite T 
and monotonic behavior of a{T = const., x) arc guaranteed. 

First we assume the transition to be continuous also at T = 0. One of the simplest 
suitable models is 



a[x 



where Tq{x Xc — 0) = 0, and a{x Xc + 0) — 0. For simplicity, all quantities are given 
in dimensionless form. 

Corresponding w{T) curves are presented in Fig. 10. Note that there are no pieces with 
w < p — 1/2 which belong to insulating samples, and that, ii w < p, dw/dT is always 
positive. 

Next we assume the transition to be discontinuous, but only at T = 0. We consider 
again a very simple qualitative model: 

l+TV2)exp{-(ro(x)/T)V2} ^<^^ 
a(x) 



-(r..)Hi;:>;r//r''' ' for (m 



where Tq{x ^ Xc — 0) =0, and a{x ^ Xc + 0) = 1. In order to illustrate how this model 
unifies a discontinuity oi a{T — const., x) at T = with continuity for T > 0, we present 
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a(T,x = const.) and a(T = const., x) in Figs. 11a and lib, respectively. The latter graph 
is based on the following additional assumptions: Xc = 0.5, Tq{x) = 30 • (xc — x)^, and 
a{x) = 1 + 3 ■ (x — Xc), where using 2 as critical exponent of Tq is motivated by p5 



however, the qualitative properties of the model considered do not depend on the numbers 
chosen. Figs. 11a and lib demonstrate two remarkable features of the model (11): (i) If 
some low-T range, e.g. the shaded area, is masked, a misinterpretation of the activated 
curves closest to the MIT as metallic ones is tempting, and one could even conclude that 
the transition is continuous, (ii) The singularity of a{0,x) at Xc evolves with decreasing T 
because a(T = const., x) becomes increasingly steep just below Xc- 

Fig. 12 shows the w{T) curves for model (11). Contrary to Fig. 10, this figure includes 
minima of w(T,x = const.) with very small values of Wmm(a;) at Tmin(x), where Tmin(x) as 
well as Wrain{x) vauish as x — Xc. The specific point is that a{0,Xc — 0) equals ct(0,Xc) for 
model (10), but not for (11). In both cases, w{0,Xc — 0) = oo, and w{0,Xc + 0) = 0; the 
limit w{T — > 0,Xc) equals 1/2 for model (10) and for model (11). 

Concerning the existence of minima in w{T,x = const.). Fig. 12 resembles far more the 
experimental figures than Fig. 10. Thus, a discontinuity of cr(0,x) is suggested, provided 
w{T, Xc) does not finally saturate (at some value far lower than theoretically expected) below 
the lowest experimentally accessible T as T decreases. However, this result was obtained 
using a very simple model. The question of its general validity is considered in the next 
subsection. 



C. Mathematical view 

The argumentation leading from the discontinuous character of the MIT, as modeled by 
Eq. (11), to the qualitative character of w{T,x) shown in Fig. 12 can easily be reversed: 
Figs. 7 and 8 suggest the assumption that w{T, Xc) ^ as T — > 0. Moreover, we only need 
the experimentally demonstrated monotonicity of a{T = const., x), that is dcr/dx > 0. We 
suppose it to be valid down to arbitrarily low T. 

Formulated in mathematical language. Figs. 7 and 8, together with Table IV, suggest 
that the samples can be discriminated according to the behavior of w{T,x) for vanishing T 
to belong to either of two groups with x < Xc and x > Xc, respectively: 

(i) For samples with "activated" conduction, classified as "insulating", there is a positive 
wi{x), such that 

w{T, x) > wi{x) for all T <T^, (12) 
where Tm is some experimentally accessible temperature. 

(ii) For "metallic" samples, there are numbers wq and q (independent of x) such that 

w{T, x)<wo- {T/T^Y for all T <T^. (13) 

We assume that we have determined (T(Tm,x) experimentally. For the "insulating" x 
region, integration of Eq. (12) leads to 

a(T,x)<a(T„,x)-(T/Tj-^(-), (14) 

so that 
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cr(0,x) = for all x < Xc- (15) 
Now we turn to the "metallic" region. For x = Xc, integration of Eq. (13) yields 

lna(T^,Xe) - \na{T,x,) <wo-[l- {T/T,,,y]/q, (16) 

so that 

cr(0,Xc) > cr(rm,Xc) •exp(-Wo/g) > 0. (17) 

That means (t(0, Xc) must be finite. (In our case, Fig. 8 leads to the estimates wq = 0.05 
and g = 1/2 for Tm = 1 K, so that exp{—wo/q} ~ 0.9.) Finally, due to the monotonicity of 
a(T = const., x), we get 

ct(0,x) > (j(Tm,Xc) ■ exp^—wo/q) for all x > Xc ■ (18) 

Hence, the classifications of the samples according to w{T —>■ 0,x) and a{T —>■ 0,x), 
respectively, are identical, and we arrive at the following conclusion: Provided our qualitative 
assumptions on the limiting behavior of w{T,Xc) and the monotonicity of cr(T = const., x) 
are justified, there must be a finite minimum metallic conductivity. 



D. Consideration of counterarguments 

We must, however, look back at our arguments to see what would be necessary for the 
transition nevertheless to be continuous, in which manner condition (13) would have to be 
weakened. Two possibilities have already been mentioned: (i) The behavior of w{T, x) might 
qualitatively change at temperatures below the lowest one accessible to us. (ii) At Xc, the 
parameter c of Eq. (4) might have a small but finite value, far smaller than the values 1/2 



and 1/3, which were predicted theoretically in Refs. P7| , p8| and used in many experimental 
studies in the literature. 

A third possibility is found when referring to Eq. (5), which, for c = 0, simplifies to 
In a{T = 0) = ln(T(Tm) — d/q. For no discontinuity to occur, it has to be possible to 
reach any (arbitrarily small) a{0,Xc + 0). The first term, lna{T^), is finite due to the 
monotonicity of cr{T^, x) and the finite a values observed in the insulating region at . 
Thus for no discontinuity to occur, d/q has to be infinite. Since d = w{T^), the exponent 
q must be infinitely small. For T independent q, this situation is equivalent to the case (ii) 
above, that is the situation of a small but finite c. The only remaining possibility is that 
of a T dependent q, vanishing as T — for x = Xc- This would imply a non-power law 
behavior of a{T,Xc)- If one would nevertheless try to fit cr oc T^, the effective exponent p 
would be T dependent, and approach as T — > 0. 

However, we know of no physical argument that suggests that the power law (which is 
derived by perturbation arguments) should weaken at low T, and we did not obtain any 
data that suggest that there is a tendency towards such a change at lowest accessible T. 
Therefore, for T = 0, although continuity across the transition is theoretically possible, 
experimental evidence appears to indicate that the transition is, in fact, discontinuous. 

Finally, the question remains whether or not this result might be an artifact arising from 
any sample inhomogeneities. We utilized all available experimental possibilities to control 
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and to minimize macroscopic inhomogeneities. Nevertheless, surely some small composition 
gradients remain. Calculating the T = conductance of samples with constant x gradients 
perpendicular and parallel to the current direction, respectively, assuming Eq. (1) to be valid, 
one obtains the following results: In the former case, the transition is "smeared out", but 
in the latter case, provided the exponent u in Eq. (1) satisfies < z/ < 1, a sharp transition 
would result. However, according to the preparation technologies used, the inhomogeneities 
should be considerably smaller parallel to the current direction than perpendicular to it. 
Thus the influence of a parallel x gradient is probably overcompensated by the perpendicular 
inhomogeneities. Moreover, the hypothesis that the discontinuous transition for T = 
might be only an artifact originating from a parallel gradient, is ruled out by the behavior 
of w{T, X = const.) for the region where a{x) ^ b{x) ■ T^. The logarithmic derivative would 
have to exhibit a plateau w(T,Xc) ~ p here according to the gradient hypothesis, but we 
did not find one. 



E. Minimum metallic conductivity 

In summary, it is very likely that a{0,x) exhibits a discontinuity in a-Sii_2,.Nii., corre- 
sponding to Mott's idea of a minimum metallic conductivity crmm tiut a(T = const. > 
0,x) should be continuous for arbitrary finite T. However, we cannot exclude the possibil- 
ity that w{T,Xc) changes its behavior, in particular that it saturates at some finite value, 
at temperatures lower than those experimentally accessible to us. This would correspond 
to the T exponent in Eq. (2) becoming much smaller than expected theoretically. Thus a 
conclusion with complete certainty is impossible. 

Our experimental data only allow us to determine bounds for cJmm. From the sample sets 
A and B we get 20 fl~^cm~^ < amm < 65 fl~^cm~^, and 25 fl~^cm~^ < cTmm < 70 fl~^cm~^, 
respectively. It is instructive to compare these data with Mott's theoretical result for cTmm 
PJT^, though electron-electron interaction and weak localization effects are neglected in 
its derivation. According to Ref. |T^, a^nm = 0.026 e^/(/ia), where a denotes the distance 



between neighboring impurity atoms, i.e., Ni atoms, at the MIT. Taking the atomic density 
of the samples to be similar that in crystalline Si, we get a ~ 5 A. This estimate fits with the 

1/3 

Mott-Edwards-Sienko criterion, • Uc ~ 0.25 with a"^ being the Bohr orbit radius of an 
isolated center and Uc the critical impurity density [[r^-|75| : Approximating by half the 
nearest neighbor distance of crystalline Ni, ^ 1.24 A^^yields just the expected product 
value. Finally, we obtain o"mm ~ 120 f2~^cm~^. This result exceeds our lower and upper 
bounds on by factors of roughly 5 and 2, respectively. Note however that the above 
consideration does not take account of the possible d character of the Ni states; moreover, 
the dimensionless parameters have some uncertainty. 

For comparison, for a-Sii-^^Cr^, the conductivity studies lead to cTmm ~ 250 Q^^cm^^ 
29| , |30[] . In that substance, annealing causes a conductivity decrease |3^, but the question 



whether or not annealing also causes a decrease of a^nm is still open. However, due to the 
differences between the phase diagrams of the systems Si-Ni and Si-Cr, see Section I, the 
following hypothesis is suggested: If, instead of the Ni-Ni distance, the distance between 
metallic NiSi2 grains were the relevant length for crmm, the relatively small value of a^m in 
a-Sii_i,Nix might result from the formation of NiSi2 crystalline regions with a diameter of 
the order of 1 to 2 nm, cf. [0,^0]. Corresponding structural studies would be interesting. 
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VII. CONCLUSIONS 



The main result of our study is the following: the T dependence of the logarithmic tem- 
perature derivative of the conductivity, w{T,x) = dlna/dlnT, for a-Sii_i.Ni^. qualitatively 
differs from the predictions of commonly accepted theory. According to this theory, cr(0, x) is 
continuous, and, on the metallic side of the MIT, (t(T, x) = a{x) + b{x) ■ T^, where p should 
equal 1/2 or 1/3. Thus w{T,Xc) = p, and, for w < p, only positive dw/dT is expected. 
However, we observed characteristic minima of w{T,x = const.), where not only T^^^ but 
also tOmin secm to tend to as x — Xc from the insulating side. Below T^in, w was observed 
to increase monotonically with decreasing T. The detailed analysis of this feature favors the 
discontinuous character of the MIT at T = 0, in the sense of Mott's original prediction of a 
finite minimum metallic conductivity. 

The question arises why a minimum metallic conductivity was found in a-Sii_a;Nia;, and 



previously in a-Sii-^^Cr^. [^^-|H3|, but not in many of the other substances, e.g. crystalline 
Si:P, investigated so far. Of course, at the present stage, one cannot exclude the possibility of 
a-Sii_^.Ni^. and Si:P belonging to different universality classes. Simultaneously, however, the 
question arises as to whether an analysis using w{T, x) of previous experiments might reveal 
similar contradictions between standard description and experimental data as we found here. 
But, for the lack of sufficiently detailed information on these investigations, this question 
has to be postponed to future studies. Such analyses should be highly interesting since 
samples with w < 1/3 and negative dw/dT at the lowest T were also found in experiments 
on crystalline Si: As and Si:(P,B) ||2T|, on amorphous Sii-ajCr^, and Sii-xMn^: ^8 



as well as on granular Al-Ge p2|. But there is also one case in the literature, where, for such 
a sample, w first increases and then clearly decreases as T vanishes: the a-Cr-SiOa, sample 
7 in Ref. [^] (see Fig. 2 of that work). Moreover, to the best of our knowledge, none of the 
experiments, which are interpreted as proof of the continuity of the MIT at T = 0, excludes 
that the "metallic" samples very close to the MIT could be insulating with a characteristic 



hopping temperature smaller than the lowest experimentally accessible T |I8|JT9[ . 

Additionally to the result concerning the character of the MIT, we consider the following 
points as important conclusions: 

- Although both the sample sets A and B, prepared with different technologies, differ signif- 
icantly concerning a{T,x), their critical Ni concentrations are close: 15.9±1 and 17±3 at%, 
respectively, related to our RBS scale. 

- A high-T mechanism, which "survives" the transition, masks the specific localization fea- 
tures to a large extent in both sample sets. It remains relevant down to temperatures 
of several hundred mK. Thus determining theoretical parameters by fitting conductivity 
formulae which account only for one mechanism is not justified for a-Sii_a;Nia, within the 
considered (T, x) region. 

- For a-Sii_a;Nia; also, it is possible to prepare metallic samples with the conductivity in- 
creasing with decreasing T. However, the infiuence of the mechanism from which this feature 
originates is weak, and it is directly visible only for certain preparation conditions. 

We would also like to stress a technological consequence of our STM investigations: When 
evaporating from two crucibles, the surface roughness of the films is large enough to lead to 
substantial concentration fiuctuations on a mesoscopic length scale, caused by fiuctuating 
incidence angles. To what extent these fluctuations are washed out by diffusion is an open 
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question. 

Finally, we turn to the question of what can be learned concerning the nature of the 
MIT in a-Sii_j;Nij. from the phenomenological results obtained above. Our findings are 
in disagreement with conclusions drawn from the treatment of the MIT in terms of one- 
parameter scaling theory 0] of Anderson localization of non-interacting electrons. But, on 
their own, they cannot explain why this description fails. For example, they do not tell 
us whether it is because electron-electron interaction had to be taken into account, nor 
whether two (or even more) scaling parameters would be needed to get an adequate result. 
Therefore, measurements of other quantities have to be consulted. XPS, UPS, and XES 
studies []50| , |5^ , as well as optical measurements have shown how the gap in the 

density of states shrinks and finally closes with increasing Ni content, and how the Fermi 
energy is shifted in this process. However, for reasons of resolution, such experiments do not 
yield information on the important immediate vicinity of the MIT since the characteristic 
energy (temperature) of the hopping conduction vanishes as x approaches Xc from below. 



A helpful hint comes from low-T hopping conduction deep in the insulating region |42,43 



- there, it has a far stronger influence on w than the high-temperature mechanism which 
survives the MIT. Since the value of the hopping exponent q in the approximation a oc 
exp{ — (Tq/T)''} is close to 1/2, and thus considerably exceeds Mott's result 1/4 for non- 
interacting electrons |80], electron-electron interaction probably is important - the question 
whether or not this value can be interpreted in terms of hopping in the Coulomb gap is still 
under controversial debate, see e.g. [pl| -|83[|. Thus, it seems likely that the MIT arises from 
interwoven localization and electron-electron interaction processes, cf. the speculations in 
Ref . M . A theoretical approach in this direction was very recently published by Chitra and 



Kotliar |8J]. These authors incorporate the long-range Coulomb interaction into dynamical 
mean-field theory, and obtain the result that the MIT should be discontinuous in two- and 
three-dimensional systems. However, if this hypothesis is true, we would be left with a new 
puzzle: Why does Mott's estimate of the minimum metallic conductivity, the derivation 
Iff^ of which neglects electron-electron interaction, still yield a reasonable value for cJmm in 



a-Sii_i:Ni 



? 

X ■ 
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TABLES 



Sample 2;/at% Ni t/nm ct(300 K)/Jl~icm~i a(4.2 K)/Jl-icm- 

1 149 159 150 15 

2 16.4 183 160 20 

3 16.7 179 245 64 

4 18.9 181 340 145 

5 22.1 162 585 330 

6 55.9 175 



TABLE I. Ni content x, film thickness t, conductivity a at room and liquid-helium tempera- 
tures for sample set A. The x values are obtained from the Ni:Si ratio measured by RBS within 
this work; their uncertainty amounts to 0.3 at% for samples 1 ~ 5, and to 0.5 at% for sample 6. 



Sample 


x/at% Ni 


t/nm 


ct(300 K)/Jl-icm-i 


ct(4.2 K)/n 


a (23) 


19.5 


121 


42 


1.0 


b (22) 


20.3 


122 


54 


2.7 


c(21) 


21.2 


124 


73 


8.0 


d (20) 


22.3 


124 


88 


18 


e(19) 


23.5 


123 


120 


30 


f (18) 


24.8 


122 


160 


60 


g(17) 


26.4 


120 


190 


82 


h (16) 


28.2 


118 


280 


160 


i(8) 


42.9 


89 


790 




j H 


25.3 









^cm ^ 



TABLE II. Ni content x, film thickness t, conductivity a at room and liquid-helium temper- 
atures for sample set B. The sample names in parentheses are the notations used in Ref. |55,56|. 



Sample j, positioned between samples f and g during film deposition, was prepared specifically for 
EDX and RBS analyses within that work. The x values are the EDX data published in either Tab. 
I of those papers. 
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Sample RBS-R RBS-F EDX-D EDX-W 

3 16?7 - 1L6 - 

3 16.7 - 12.0 

6 55.9 - 51.7 

f - - 14.0 24.8 

h - - 20.3 28.2 

i - - 37.8 42.9 

j - 20.0 - 25.3 

TABLE III. Comparison of Ni contents as determined by RBS in Rossendorf (RBS-R) and 
Faure (RBS-F), and by EDX in Dresden (EDX-D) and Witwatersrand (EDX-W). The RBS-R and 
EDX-D data were obtained within this work, whereas the RBS-F and EDX-W results are taken 



from Ref. |^^, - unfortunately, Ref. |55| contained a mistake in evaluating x from the Si:Ni ratio 
obtained by RBS-F, which is corrected here. All Ni contents were obtained from the Si:Ni ratios, 
ignoring impurity contamination. The random error of the EDX measurement within this study 
does not exceed 0.5 at% Ni. The first two lines serve also as reproducibility check of EDX-D; these 
measurements were performed fully independently at two different days, when different pieces of 
the film were investigated. 



q 


c 
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w{T 0) 


a{T - 0) 
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— CXD 
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insulator 
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ideal metal 
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finite 


real metal 
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> 





insulator 
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arbitrary 
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ideal metal 
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arbitrary 
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finite 


real metal 


> 


> 


arbitrary 


> 





insulator 



TABLE IV. Behavior of (t(T — > 0), and character of the conduction process in dependence on 
the parameters in Eq. (4), modelling w{T 0). Here, "ideal metal" stands for a{T ^ 0) = oo (no 
impurity scattering). For simplicity, if q = 0, we assume d = 0; moreover, we consider d = only 
for q = 0. 
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FIGURES 



FIG. 1. Experimental parameter plane: The insulating region, where only activated conduc- 
tion occurs, is marked by I, the metallic region by M, and the lowest experimentally accessible 
temperature by Tica- Measuring a{T,x = const.) means to obtain data points (•, x) along vertical 
lines in this graph. The characteristic hopping temperature, Thop(a;), is marked by a full line. Only 
for T < Thop, cr depends exponentially on T. For T > Thop, comparably flat, non-exponential cr(T) 
dependences are expected. Thus, for •, depending on T, both exponential and non-exponential 
behavior are observed. However, for x, only non-exponential a(T,x = const.) are found, although 
this sample belongs to the insulating region, too. The latter problem occurs in the whole interval 

[x*,Xc). 

FIG. 2. Glow discharge optical emission spectroscopy (GDOES) depth profiles for the elements 
B, C, Ni, and Si, i.e., intensity I versus sputtering time t, for two samples prepared by technologies 
A and B, respectively, (a) presents the analysis of sample 5, whereas (b) shows data for a sample 
prepared specifically for this analysis together with samples a - i; due to the simultaneous prepa- 
ration, (b) is representative for the whole sample set B. We made use of the following lines: B - 
208.9 nm, C - 156.1 nm, Ni - 349.2 nm, and Si - 288.1 nm. The intensity scales were adjusted to 
reach a high resolution, thus the units differ from element to element. In (a), the sudden increase 
of the B signal indicates the moment when the bottom of the sputtering crater reaches the glass 
substrate. 

FIG. 3. STM image of the surface topography of sample 3, measured with a current of 0.5 nA 
and a bias of 8 V. The figure was obtained using a slope correction and a slight smoothing of the 
row data. 

FIG. 4. Three STM cross sections of the surface of sample 3 obtained with different resolutions 
(current = 0.3 nA, bias = 8 V). The data presented are obtained from the row data performing 
only a slope correction. In all cases, the height was determined at 128 points along a line segment, 
slightly longer than represented. 

FIG. 5. Overview of a{T,x) in double logarithmic representations for sets A (+) and B (x) 
in graphs (a) and (b), respectively. For the characterization of the samples see Tables I and II. 
To illustrate common as well as differing features, both graphs include a sample of the other set, 
for the sake of distinctness marked by •. Moreover, (a) includes data for two sputtered a-Sii^^^Nij; 
films (a) from Ref. |^^. These samples, here referred to as a and (3, have a Ni content of 20.7, 
and 27.1 at%, respectively, as determined by an EDX analysis |50[| . 

FIG. 6. a{T,x) versus T^/^ representation of the samples closest to the MIT for sets A (+) 
and B (x) in graphs (a) and (b), respectively. Both graphs include one sample of the other set 
(•). The dashed lines give the extrapolations according to Eq. (2) with p = 1/2, obtained from the 
region 2 - 30 K. 
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FIG. 7. Overview of the T and x dependence of the logarithmic derivative, w = dlncr/dlnT 
for sets A (+) and B (x) in graphs (a) and (b), respectively. Both graphs include one curve of the 
other set (•). The temperature scale T^/^ is chosen in order to have reasonable resolution also for 
the low-T part. 

FIG. 8. Low-T / low-w part of Fig. 7. The data presented were obtained from measurements 
in three different cryostats, see Section IV, where several months passed between the investigations. 
Moreover, in some cases, contacts had to be renewed between the measurements. 

FIG. 9. Comparison of w{T) for the a-Sii_j;Nia; sample 4 (•) with the w{T) for three a-Sii-^Crj; 
samples (*), here referred to as 7, 6, and e. The a-Sii-xCr^; curves were obtained from cr{T) data, 
published in Figs. 1 and 7 of Ref. |^^. Samples 7, 5, and e, have a Cr content of 11, 14, and 
19 at%, respectively, obtained by EDX |3^. Calculating w, we used /c = 4 for samples 7 and e, 
and k = 2 for 5. 

FIG. 10. Qualitative behavior of w{T,x) for a continuous transition according to Eq. 

(10) . Dashed lines: insulating with Tq{x)^^'^ = 0.03, 0.1, 0.3, and 1; full lines: metallic with 
a{x) = 0.03, 0.1, 0.3, and 1; dashed-dotted line: separatrix between insulating and metallic re- 
gions. 

FIG. 11. Characteristic a{T,x) features of a discontinuous transition according to Eq. (11): 
T and x dependences of a are shown in graphs (a) and (b), respectively. Dashed and full lines 
represent insulating and metallic regions, respectively; the dashed-dotted line in (a) denotes the 
separatrix between these regions. The lines shown in (a) were obtained for the parameter values 
To{xy/'^ = 0.05, 0.15, 0.5, 1.5, and a{x) - 1 = 0.2, 0.6, 1.0, 1.4. Shaded area in (a): example of 
low-T range being inaccessible in a real experiment. 

FIG. 12. Qualitative behavior of w{T,x) for a discontinuous transition according to Eq. 

(11) . Dashed lines: insulating with To(x)^/^ = 0.02, 0.06, 0.2, and 0.6; full lines: metallic with 
a{x) — 1 = 0.2, 0.6, 2, and 6; dashed-dotted line: separatrix between insulating and metallic 
regions. 
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Fig. 2a 




Fig. 2b 



This figure "fig_03.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/98 1 2 106v3 
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